Emergent Exclusion Statistics of Fibonacci Anyons in 2D Topological Phases 
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We demonstrate how the generalized Pauli exclusion principle emerges for quasiparticle excitations 
in 2d topological phases. As an example, we examine the Levin- Wen model with the Fibonacci 
data (specified in the text), and construct the number operator for fluxons living on plaquettes. 
By numerically counting the many-body states with fluxon number fixed, the matrix of exclusion 
statistics parameters is identified and is shown to depend on the spatial topology (sphere or torus) 
of the system. Our work reveals the structure of the (many-body) Hilbert space and some general 
features of thermodynamics for quasiparticle excitations in topological matter. 
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1. Introduction 

By now it is well-known that (quasi-)particles in 
strongly entangled many-body systems may exhibit ex- 
otic quantum statistics (see, e.g., ref. [l| for a review), 
other than the usual Bose-Einstein and Fermi-Dirac 
ones. In addition to the anyonic or braiding exchange 
statistics 2 -^ in two dimensional systems, statistical weight 
of many-body quantum states may also obey new com- 
binatoric counting rules^^, in which the number of avail- 
able single-particle states, when adding one more quasi- 
particle into the system, linearly depends on the number 
of existing quasi-particlcs. A typical new feature of the 
generalized Pauli exclusion principle is mutual exclusion 
between different species, resulting in a matrix of statis- 
tical parameters^, as well as unusual thermodynamics for 
ideal gases with only statistical interactions^. 

More precisely, following ref. [H, in the case with only 
one species of quasi-particles, the number of A^-particlc 
states is assumed to be given by the binomial coefficient: 
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with G c ff = G — a(N — 1) being the number of available 
single-particle states, while G is the number of single- 
particle states when N = 1. Then a = corresponds to 
bosons and a = 1 fcrmions; other values of a gives rise to 
exotic exclusion statistics. Similarly, in the multi-species 
case, the number of many-particle states is assumed to 
be given by (a, b = 1, . . . , m labeling species) 
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Here coefficients a a b form the (mutual) statistics matrix. 

It has been shown£ that the Thermodynamic Ansatai 
for one-dimensional solvable many-particle models is ac- 
tually a special case of the exotic exclusion statistics. 
(See also refsHH.) It has been also numerically verified 
that quasi-particle excitations in the fractional quantum 



Hall (FQH)systems indeed obey^ eq. (TJ, or eq. © 
allowing mutual exclusion between different species^. 
Moreover either the Haldane or Jain hierarchy in the 
FQH effect can be theoretically understood from the ex- 
clusion statistics of quasiparticle s 1 ' 1 2 . 

Recently there has been revived interest in the study 
of quasiparticle statistics in 2d topological states of mat- 
ter (including FQH systems), because of the possibility 
of using their braiding to do (fault tolerant) topologi- 
cal quantum computation (TQC) 13 ' 14 . In order to know 
better about the error of TQC at finite temperature, it is 
needed to understand better how exclusion statistics of 
quasi-particles emerges in 2d topological matter, which 
governs the thermodynamics of the system. 

In this paper, we carry out the many-body state count- 
ing in an exactly solvable discrete model, i.e. the Levin- 
Wen modet^ (with a special set of data) , that describes 
a 2d topological quantum fluicU& of Fibonacci anyons^, 
with doubled Fibonacci anyons as fluxon excitations liv- 
ing on plaquettes. The Fibonacci anyons are the simplest 
non-abelian anyons. They occur as quasiparticles in the 
k = 3 Rcad-Rezayi stated in a FQH state with filling 
fraction v = 4?-, and can be used for universal topologi- 
cal quantum computation^. (Recently, it is proposed^ 
that the physics of interacting Fibonacci anyons may be 
studied in a Rydbcrg lattice gas.) 

In this Letter, we first construct the number opera- 
tor for fluxons in the model, which helps us identify 
the states with localized excitations. Then we numeri- 
cally count the (many-body) states with fluxon-numbcr 
N fixed, from = 1 up to N — 7, for the system on 
a sphere and torus, respectively. The results exhibit a 
pattern closely related to the Fibonacci numbers, which 
in turn is put in the form of eq. thus determining 
a topology-dependent statistical parameter matrix. Our 
work reveals that exotic exclusion emerges among quasi- 
particlcs due to interplay between various "hidden" de- 
grees of freedom in addition to fluxon locations. Finally, 
we briefly discuss the thermodynamics of the system. 

2. The model 
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We consider a discrete model for a "spin" system on 
a trivalent graph on a closed surface, e.g. a sphere or 
torus. We adopt a simplified formulation of the Levin- 
Wen models, with the Fibonacci data (e.g., see ref. [13 ) 
given below. Each link is assigned a "spin" -type j = 
0, 1, and configurations of the labels on all links form an 
orthonormal basis in the Hilbert space. (It is conjectured 
that the Levin- Wen model describes a class of doubled 
(time-reversal invariant) topological phases^.) 

The Hamiltonian of the model is 



s=0,l 



where the summation runs over vertices v and plaqucttcs 
p, with d a = Mi =<(>= (v / 5+ l)/2, and D = 1 + 2 . U 
and e are positive constants. 

The operator Q v defined at vertex v is 



(4) 



(Only the relevant part of the graph is shown; the rest 





(5) 



remains the same.) Here 8^ = Sjki = Sjik is given by 

^ooo = ooii = ^111 = 1; #ooi = 

(called the Fibonacci fusion ruleii). 
The operator Bp defined at plaquette p is 



Bi 



= Vjl v h v i3 v r v» v r (:'■! '/■ C;' ; / '/ CP 9 ?,' i 2 , (6) 

01 J2 03 J 1 ]. 2 3 3 sj'sti S3[j' 2 sj' 2 j' 3 ' V / 

where Vj = \fdj. A similar rule applies when the plaque- 
tte p is a quadrangle, a pentagon, or a hexagon etc. Note 
that the matrix is nondiagonal only on the labels of the 
boundary links (i.e., j\, j2, and J3 on the above graph). 
The non-zero 6j-symbols G are given by 
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^000 
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with the (tetrahedral) symmetry: 
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One can check that they satisfy the conditions: 
pentagon id: £ d^G^A™ = G%G 

71=0,1 

orthogonality: ]T d n G^G l p Tn = 5 fS mlq S kip . (9) 

n=0,l 1 

A notable property of the model is that Q v and B p are 
mutually commuting projection operators: [Q v ,B p ] = 0, 




FIG. 1: A configuration in the subspace with only fluxons 
allowed, with solid lines for j = 1 and dotted lines j = 0. 
Open strings (with only single solid line at some vertex) are 
forbidden by 5ooi = 0. 
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and BpBp> 



dpp' Bp . 



Thus the Hamil- 
tonian is exactly solvable. The energy eigenstates are 
the simultaneous eigenvectors of these projections Q v 
and Bp. The ground states are those |<3?) that satisfy 
Q v \&) = |3>) = B „ | $ ), for all v and p. Using the method 
developed in ref. [2l|, one can compute ground state de- 
generacy: GSD = 1 on sphere and GSD = 4 on torus. 

The quasiparticle excitations are the states with zero 
eigenvalue of Q V ' for some v' and/or of B p i for some p' . 
In this paper we restrict ourselves to study the so-called 
fluxons, satisfying Q v = 1 for all v, B p i = for a specified 

set of p' (where fluxons live), and B p = 1 for all other p. 
(See Fig. H) 

3. Number operator of fluxons 

A crucial property of the model is that Bp defined 
above form an abelian algebra: 



b p B p = 



(10) 



which can be verified using conditions (JSJ) and ©• From 
the generators we construct the operators (i = 0, 1): 



where S is the modular matrix given by 



S = 



1 



1 



1 



(11) 



(12) 



with Sij for fixed i being a one-dimensional representa- 
tion of the algebra (fT0|) . One can check that n p (i = 0, 1) 
form a complete set of orthonormal projections: 



nln 1 ! = 6 ik n 



nl = l. 



(13) 



There is a fluxon at p in a state |^), if np\^) 
A ground state |$) contains no fluxon because n 



|*). 
B p . 

Hence the model has only one type of fluxons. No state 
contains two fluxons at any p. However, we will show 
that the fluxon in this model is not a fermion; rather its 
state counting needs exotic exclusion statistics. 

4. Exclusion Statistics On a sphere 

Let us count the iV-fluxon states in the model with P 
plaquettes on a sphere. Pick up a set of N fixed pla- 
quettes and denote it by C — {p\,p2, ■■■,Pn} {N < P). 
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The states with exactly N fluxons occupying the selected 
plaquettes are those satisfying 

^l*) = M*). forp'^C. (14) 

Thus ^n p ec Ilp'^c "p') * s ^ nc projector onto the 
subspace of such states. Tracing this projection computes 
the total number of the A^-fluxon states in the configu- 
ration C: 

«;p,jv,c=tr(n^n^)- ( is ) 

pec P '£c 

We numerically compute eq. (|15p on random graphs 
on sphere with P(> 7) plaquettes, with the stable result 
presented in Table [T] 



TABLE I: State Counting on Sphere 



Fluxon number N 





1 


2 3 


4 


5 


6 


7 


State Counting wp t N,c 


1 





1 1 


4 


9 


25 


64 



The pattern of the iV-dependence is obvious: 

Wp,n,C = Fn-v ( 16 ) 

where F n is the Fibonacci number that satisfies the re- 
currence relation: F n = -F n _i + -Fri-2 with F\ = Fi = 1. 
Both numerically and analytically we have checked that 
eq. (|16|) is independent of the graph, of the total number 
P of plaquettes, as well as the locations of the N fluxons. 
The appearance of the squared in eq. (|16p is consistent 
with the conjecture that the LW model describes a dou- 
bled topological phases^. 

Summing over configurations C (i.e. over possible dis- 
tributions of N plaquettes in a fixed graph) , we get the 
total number of iV-fluxon states: 

iySr=E tt w = (Jk-i- (17) 

The first factor counts the ways to distribute A" fluxons 
over P plaquettes. The second factor counts the config- 
urations of the link degrees of freedom, which are not 
unique, given N and C. The independence of wp t N,c on 
P and C implies the degeneracy of the excited states is 
topological in the sense that it does not depend on the 
detailed structure of the underlying graph, and not on 
the relative positions between the fluxons as well. The 
origin of this property lies in the topological symmetry 
of the model under mutations of the underlying graphSl. 

To find the exclusion statistics, we rewrite p7|) : 

v 7 JVi,iV2=0 v ' v y 

(18) 



where [x] is the greatest integer less than or equal to x. 

Now eq. @ is of the form of eq. (fT5|) . by introduc- 
ing two additional pseudo-species a = 1,2, which do 
not contribute to the total energy but are helpful for 
state-counting. This is similar to what was suggested 
for state counting in some conformal field theories^. In- 
cluding the original fluxon species labeled by a = 0, from 
cq. (|18[) we read the exclusion statistics parameters a a t 
(a, b, = 0,1, 2): 



-12 0. (19) 



The diagonal a aa is the self-exclusion statistics for 
species a. The aoo = 1 implies the hard-core boson be- 
havior. This can be understood with cq. fj 13|) . 

The pseudo-species provides a way to count config- 
urations, in the presence of fluxons, of link degrees of 
freedom, which are not uniquely determined by the con- 
straints (fT4|) . The value an = OL21 = 2 implies that 
one pseudo-particle makes two single-particle states (or 
"seats") unavailable to an additional pseudo-particle. 
The negative mutual statistics «20 = Q^30 = ~ 1 tells us 
that each fluxon present creates one vacant "seat" for 
each pseudo-species. So the maximum particle number 
of each pseudo-species is naturally [(N — l)/2]. These 
results help us to understand the structure of the (many- 
body) Hilbert space for excited states of the system, and 
to derive analytically the state counting formula (|18p. 

6. Exclusion statistics on a torus 

We proceed to consider the model on a torus. The 
ground state degeneracy^ is 4 . Thus the system ex- 
hibits the global topological degrees of freedom, and we 
can study their effects on excited states by counting the 
pseudo-particle states. 

Pick up N plaquettes (N < P). The number of states 
with fluxons on htesc plaquettes is computed numeri- 
cally as in Table [HI 



TABLE II: State Counting on Torus 



Fluxon number N 





1 2 


3 


4 


5 6 


State Counting 


2 2 


1 3 2 


4 2 


7 2 


ll 2 18 2 



The pattern of its dependence on A^ is 

w£T=(£W> (20) 

with L n the Lucas number, a modified version of the Fi- 
bonacci number, satisfying the recurrence relation L n = 
L n -i + L„_ 2 with Li = 1,L 2 = 3. 
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We rewrite (l20l) in terms of binomial coefficients: 



totically 



v 7 iV 1; iV2=0,l x 7 v ^ 



[i(W-2)] 

E 

^3,^4=0 



iV — 27Vi — Y 3 \ /AT - 2Y 2 - Y 4 
Y 3 J 1 N 4 



(21) 



and get the exclusion statistics parameters a a b (o>, b = 
0,1,2,3,4): 



a 



( 1 o o\ 
10 
10 
-f 2 2 
V -1 2 2 y 



(22) 



where we denote by a = the ffuxon species. 

Eq. (f2"Tj) shows that one needs to introduce four 
pseudo-species a = 1, 2, 3, 4. The pseudo-species a = 1, 2 
are interpreted as the topological degrees of freedom 
on the torus, for the following reasons. The allowed 
"particle number" Ni,N2 — 0,1 of these pseudo-species 
are independent of the number AT of fiuxons. Particu- 
larly when there is no fluxon present, the configurations 
N\,N2 = 0,1 characterize the four-degenerate ground 
states. Then the pseudo-species a = 3,4 provide a way 
to count the configurations of link degrees of freedom 
given a ground state and fluxon number. 

The state counting of excitations on a torus is shown 
different from that on a sphere. Indeed the mutual statis- 
tics parameters a 3 i = a±2 = 2 imply that the number of 
configurations of link degrees of freedom a = 3 (a = 4) 
are affected by the topological degrees of freedom a = 1 
(a = 2), respectively. On the other hand, the topolog- 
ical degrees of freedom are not affected by the fiuxons 
present and the link degrees of freedom. So the degener- 
ate ground states can be used to label the sectors of ex- 
citations. We note that in the sector with N% = N% = 1, 
the state counting for fiuxons is exactly the same as that 
on sphere. 

7. Statistical Thermodynamics 

Now we assume that only fiuxons can be thermally ex- 
cited; this is the case when U » e,kT in eq. In 
the thermodynamic limit, the Hilbert space dimension of 
Y-fiuxon states (occupying N fixed plaqucttcs) is asymp- 



on torus: 



lim L A 



N-yoo 



N 



on sphere: lim F^-_ x - <t> 2N ~ 2 /5, 

i2N 



(23) 



(cf) 2 is called the quantum dimension of the fluxon.) On 
a torus, for example, the canonical partition function is 



E 

N=0 



L 2 N e 



-Ne/kT 



b 2 e -e/kT + l] 



(24) 



It can be interpreted as the grand canonical partition 
function of the many-fluxon system, which behaves like a 
fermionic system with a temperature-independent fugac- 
ity z given by the quantum dimension: 



(25) 



The fugacity z counts the effective number of states per 
fluxon located at a plaquettc. Note that z is irrational 
rather than integer. This is a manifestation that the 
many-fluxon states are highly entangled ones with long- 
range entanglement. They are superpositions of highly 
constrained j-configurations on the links, obviously not 
of the form of a direct product of localized fluxon states. 

The statistical distribution of the average occupation 
number of fiuxons is obtained from eq. (|24|) : 



(n) = (N)/P = 



1 



1 



(26) 



Many useful thermodynamic observables are then com- 
putable. Though the model is very simple, we believe 
that the features revealed in this paper should be quite 
general for emergent exotic exclusion statistics and ther- 
modynamics for quasiparticle excitations in a wide class 
of 2d topological phases. 
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Note added. After this paper was posted on arXiv, 
we learn that a different state-counting formula, without 
proof and connection to exclusion statistics, appeared in 
M. D. Schulz, S . Dusuel, K. P. Schmidt, and J. Vidal, 
larXiv:1212.4109l 
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